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Abstract 

The theory of integral quadratic constraints (IQCs) allows verification of stability 
and gain-bound properties of systems containing nonlinear or uncertain elements. Gain 
bounds often imply exponential stability, but it can be challenging to compute useful 
numerical bounds on the exponential decay rate. In this work, we present a modifica¬ 
tion of the classical IQC results of Megretski and Rantzer [13] that leads to a tractable 
computational procedure for finding exponential rate certificates. We demonstrate the 
effectiveness of our method via a numerical example. 


1 Introduction 

In robust control problems, we seek absolute performance guarantees about a system in 
the presence of bounded uncertainty. Examples of such results include the small gain the¬ 
orem & passivity theory [20], dissipativity theory [19], and integral quadratic constraints 
(IQCs) [13]. 

In this paper, we present a modification of IQC theory, the most general of the afore¬ 
mentioned tools, that allows one to certify exponential stability rather than just bounded- 
input bounded-output (BIBO) stability. Moreover, we can compute numerical bounds on 
the exponential decay rate of the state. 

Even when BIBO stable systems are exponentially stable, the estimates of the expo¬ 
nential decay rates provided by standard IQC theory are typically very conservative. We 
will show that this conservatism can be greatly reduced if we directly certify exponential 
stability and use the method presented herein to compute the associated decay rate. 

Our modified IQC analysis was successfully applied in [12] to analyze convergence prop¬ 
erties of commonly-used optimization algorithms such as the gradient descent method. 
These algorithms converge at an exponential rate when applied to strongly convex func¬ 
tions, and the modified IQC analysis automatically produces very tight bounds on the 
convergence rates. 

Another potential application is in time-critical applications such as embedded model 
predictive control, where it is vital to have robust guarantees that the desired error bounds 
will be achieved in the allotted time without overflow errors and in spite of fixed-point 
arithmetic. See [8] and references therein. 


A special case. As previously noted, exponential stability certificates are often con¬ 
servative when they are derived from L 2 gain bounds. However, it is well known that 
exponential stability can be proven directly in some special cases. To illustrate this 
fact, consider a discrete linear time-invariant (LTI) plant G with state-space realization 
(A, B,C, D). Suppose G is connected in feedback with a passive nonlinearity A. A suf¬ 
ficient condition for BIBO stability is that there exists a positive definite matrix P >- 0 
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and a scalar A > 0 satisfying the linear matrix inequality (LMI) 
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If we define V(x) := x T Px, then (1) implies that V decreases along trajectories: 
V(xk+i) < V(xk) for all k. BIBO stability then follows from positivity and bound¬ 
edness of V. But observe that when (1) holds, we may replace the right-hand side by 
—eP for some e > 0 sufficiently small. We then conclude that V(xk+ 1 ) < (1 — e)V(xk) 
for all k and exponential stability follows. We can then maximize e subject to feasibility 
of (1) to further improve the rate bound. 

Unfortunately, the simple trick shown above does not work in the general IQC setting 
due to the different role played by P in the associated LMI. The LMI used in IQC theory 
comes from the Kalman-Yakubovich-Popov (KYP) lemma and although it is structurally 
similar to (1), P is not positive definite in general and V may not decrease along trajec¬ 
tories. 

Our key insight is that with a suitable modification to both the LMI and the IQC 
definition, we obtain a condition that can certify exponential stability. 

The paper is organized as follows. We cover some related work in the remainder of 
the introduction, we explain our notation and some basic results in Section 2, we develop 
and present our main result in Section 3, and we discuss computational considerations 
in Section 4. Finally, we present an illustrative example demonstrating the usefulness of 
our result in Section 5, and we make some concluding remarks in Section 6. 


Related work. It is noted in [13, 15] that BIBO stability often implies exponential 
stability. In particular, we get exponential stability if the nonlinearity satisfies an addi¬ 
tional fading memory property. The proof of this result is chiefly concerned with showing 
existence of an exponential decay rate. Although the proof constructs an exponential 
rate, the construction is based on the assumed L 2 gain of the linear map, and thus can 
be very conservative. 

Other proofs of exponential stability have appeared in the literature for specific classes 
of nonlinearities. Some examples include [4, 11], which treat sector-bounded nonlineari¬ 
ties, and [9], which treats nonlinearities satisfying a Popov IQC. These works exploit LMI 
modifications akin to the one shown with (1) earlier in this section. 

The sequel is inspired by the recent paper [12], which presents an approach for proving 
the robust exponential stability of optimization algorithms. The approach of [12] uses a 
time-domain formulation of IQCs modified to handle exponential stability. In contrast, 
the present work develops the aforementioned exponential stability modification entirely 
in the frequency domain and clarifies its connection to the seminal IQC results in [13]. 

2 Notation and preliminaries 

We adopt a setup analogous to the one used in [13], with the exception that we will work 
in discrete time rather than continuous time. The conjugate transpose of a vector v £ C n 
is denoted v*. The unit circle in the complex plane is denoted T := {z € C | \z\ = 1} 
The ^-transform of a time-domain signal x := (xo, x\,. .. ) is denoted x(z) and defined as 

x(z) := YX=o x k z ~ k - 

A Hermitian positive definite (semidehnite) matrix M is denoted M >- 0 ( M Y 0). 
Function composition is denoted (go f)(x) := g(f(x)). A sequence u = (uq,Ui, ...) is 
said to be in £2 if SfcLo l Ufc l 2 < 00 . A sequence Uk is said to be in t p 2 for some p € (0,1) 
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if the sequence ( p~ k Uk ) is in l 2 , he. YlkLo P~ 2k \ u k\ 2 < oo. Note that t p 2 C l 2 . Let 
IZTCffifi 11 be the set of to x n matrices whose elements are proper rational functions with 
real coefficients analytic on the closed unit disk. 

Consider the standard setup of Fig. 1. The block G contains the known LTI part of 
the system while A contains the part that is uncertain, unknown, nonlinear, or otherwise 
troublesome. The interconnection is said to be well-posed if the map (n, w) K > (e, /) has a 
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Figure 1: Linear time-invariant system G in feedback with a nonlinearity A. 

causal inverse. The interconnection is said to be BIBO stable if, in addition, there exists 
some 7 > 0 such that when G is initialized with zero state, 

lbl| 2 + lkl| 2 <7(l|e|| 2 +ll/ll 2 ) 

for all square-summable inputs f and e, and where ||-|| denotes the i 2 norm. Finally, the 
interconnection is exponentially stable if there exists some p £ (0,1) and c > 0 such that 
if / = 0 and e = 0, the state Xk of G will decay exponentially with rate p. That is, 

\\xk\\ < cp k \\x 0 \\ for all k. 


We now present the classical IQC definition and stability result, which will be modified 
in the sequel to guarantee exponential convergence. These results are the discrete-time 
analog of the main IQC results of Megretski and Rantzer [13]. 

Definition 1 (IQC). Signals y £ t 2 and u £ £ 2 with associated z-transforms y(z) and 
u(z) satisfy the IQC defined by a Hermitian complex-valued function II if 


y{z) 

H z ). 


n(z) 


y ( z ) 

w(z) 


dz > 0 . 


( 2 ) 


A bounded operator A satisfies the IQC defined by II if (2) holds for all y £ in with 
u = A (y). We also define IQC(II(z)) to be the set of all A that satisfy the IQC defined 
by n. 

Theorem 2 (Stability result). Let G(z ) £ IZ'Hff > xn and let A be a bounded causal oper¬ 
ator. Suppose that: 


i) for every r £ [0, 1] , the interconnection of G and r A is well-posed. 

ii) for every r € [0, 1], we have rA £ IQC(n(z)). 

Hi) there exists e > 0 such that 


\G(z)V 

I 


n(z) 


G{z) 

I 


A — el, 


Vz £ T. 


Then, the feedback interconnection of G and A is stable. 
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3 Frequency-domain condition 

In this section, we augment Definition 1 and the classical result of Theorem 2 to derive a 
frequency-domain condition that certifies exponential stability. 

Definition 3. The operators p + , p_ are defined as the time-domain, time-dependent 
multipliers p k ,p~ k , respectively, where p G (0,1) is a defined constant. 

Remark 4. The operator p- o (G(z) o p + ) is equivalent to the operator G(pz). This 
follows from the fact that, for any constant a > 0 and signal Uk, the z-transform of 
a~ k Uk is given by u(az). See Fig. 2 for an illustration. 


u(z) u(p 1 z) G(z)u(p x z) G(pz)u(z) 



Figure 2: Illustration of Remark 4. 


In order to show exponential stability of the system in Fig. 1, we will relate it to BIBO 
stability of the modified system shown in Fig. 3. This equivalence is closely related to 
the theory of stability multipliers [16]. 
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Figure 3: Modified feedback diagram with additional multipliers and inputs. For appro¬ 
priately chosen e and / and with zero initial condition, we show how this diagram is 
equivalent to that of Fig. 1. 


Proposition 5. Suppose G(z) has a minimal realization (A, B,C, D). If the intercon¬ 
nection in Fig. 3 is stable with zero initial condition, then the interconnection in Fig. 1 
with initial state Xq is exponentially stable. 


Proof. Intuitively, if v and w are small in the BIBO sense compared to e and /, then y 
must be even smaller. A complete proof is included in the appendix. ■ 


Ideally, we would like to find a suitable redefinition of the IQC for this transformed 
system shown in Fig. 3. To this end, we introduce the concept of the p-IQC. 

Definition 6 (p-IQC). Signals y € t p 2 and u £ with associated z-transforms y(z) and 
u(z) satisfy the p-IQC defined by a Hermitian complex-valued function II if 



y{pz ) 

u{pz) 


dz >0 . 


( 3 ) 
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A bounded operator A satisfies the p-IQC defined by II if (3) holds for all y £ w ith 
u = A (y). We also define IQC(n(,z), p) to be the set of all A that satisfy the p-IQC 
defined by II. 

Note that the concept of a p-IQC generalizes that of a regular IQC. Indeed, we have 
IQC(II(z),l) = IQC(II( 2 :)). The restriction of u G if) an d V G i) corresponds to the 
restriction of u £ ti and y £ 1 2 in the classical definition of IQC [13]. Now equipped with 
p-IQCs, we can relate A' in Fig. 3 to A in Fig. 1. 

Proposition 7. Let A be a nonlinearity, and let II be a Hermitian complex-valued func¬ 
tion. As in Fig. 3, define A' := p_ o (A o p + ). Then the following statements are 
equivalent. 

(i) A £ IQC(n( 2 ),p) 

(it) A' £ IQC(n(pz)) 

Proof. We define the discrete Fourier transform of the input and output of A as y(z) and 
u(z), respectively. Then, from the definition of p+ and p_, we have that w(z) = u(pz) 
and v(z) = il(pz). Substituting into the IQC definition (2), we obtain (3) as required. ■ 

Proposition 7 is illustrated in Fig. 4. 
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Figure 4: Illustration of Proposition 7. 


We now state our main result, an exponential stability theorem analogous to the clas¬ 
sical result in Theorem 2. 

Theorem 8 (Exponential stability). Fix p £ (0,1). Let G(pz) £ IZTCff^ 71 and let A be a 
bounded causal operator. Suppose that: 


i) for every r £ [0,1] , the interconnection of G and r A is well-posed. 

ii) for every r £ [0,1], we have r A £ IQC(n(,z), p). 

Hi) there exists e > 0 such that 


G(pz) 

I 


n (pz) 


'G(pz) 

I 


A —eT, 


VzeT. 


(4) 


Then, the interconnection of G and A shown in Fig. 1 is exponentially stable with rate p. 

Proof. Roughly, we apply Theorem 2 to the interconnection in Fig. 3 with operators G' 
and A' and IQC II(pz). 

(a) Since Fig. 1 and Fig. 3 have the same interconnection structure, well-posedness is 
equi^valent. 
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(b) Due to the equivalence of IQCs in Proposition 7, 


t A G IQC(n( 2 ),p) 

p_ O ((rA) o p + ) g IQC(n(pz)) 

^==> r(p_ o (A o p + )) G IQC(n(p 2 )) 
rA' g IQC(n(pz)) . 

(c) This is condition iii) of Theorem 2 using G' and A'. 

Thus, these three conditions ensure BIBO stability of the system in Fig. 3. We then apply 
Proposition 5 to arrive at exponential stability of Fig. 1. ■ 


4 Computation 


As in the classical IQC setting, to guarantee stability, the frequency-domain inequality 
(FDI) (4) must be verified for every w G [0,27r). However, if the IQC in question exhibits 
a factorizaton , then the discrete-time KYP Lemma can be applied to convert the infinite¬ 
dimensional FDI to a finite-dimensional LMI. We now review these results. 

Definition 9. We say n has a factorization (4/, M) if 

n(» = $(z)*M$(z) , 

where 4/ is a stable linear time-invariant system, M is a constant Hermitian matrix, and 
^(z)* denotes the conjugate transpose of’i'(z). 

Remark 10. IfU(z) has a factorization (^, M ) and ^t(pz) is stable, then (3) is equivalent 
to 

oo 

p~ 2k zjMzk > 0 , where z := *F 

k =o 

This follows immediately from ParsevaTs theorem. 



The KYP lemma, stated below, is attributed to Kalman, Yakubovich, and Popov. A 
simple proof and further references can be found in [14]. 

Lemma 11 (Discrete-time KYP Lemma). Given matrices A,B and a Hermitian matrix 
M, and assuming A has no eigenvalues on the unit circle, the FDI 

\zI-A)-'B\* 

I 


(zl - A)~ 1 B 
I 


-< 0 


holds for all z G T if and only if there exists a solution P = P T and A > 0 to the LMI 

+ XM <0. 


A B 

T 

P o ‘ 


A B 

I 0 


0 -P 


I 0 


Corollary 12. Suppose the realization of G is given by (A, B,C, D) and assume n has 
a factorization where the realization of ’F is given by 


Aq, 

P'1'1 

B \Jr 2 

cv 

D^! 1 
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Then (4) is equivalent to the existence of P = P T and A > 0 such that 


' A t pa - p 2 p 

a t pb 
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C T 

b t pA 

b t pb 

D T _ 


M [C D\ < 0 


( 5 ) 


where (A, B , C, D) are given hy 
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Proof. A similar result is proven in [18], which we repeat here for completeness. 
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where * denotes the repeated part of the quadratic form. Similarly, we have 
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If p~ l A has no eigenvalues on the unit circle, we may then invoke Lemma 11 (applied to 
p~ 1 A 1 p~ 1 B 1 and the appropriate M term) and multiply through by p 2 to show that (4) 
is equivalent to the existence of P = P T and A > 0 such that (5) holds, as required. ■ 


With the advent of fast interior-point methods to solve LMIs, the feasibility of the LMI 
(5) can be quickly ascertained for any fixed p 2 . Since the size of the LMI is on the order 
of the size of the system G and the IQC II, most practical linear systems lead to relatively 
small LMIs. 

Finding the best upper bound amounts to minimizing p 2 subject to (5) being feasible. 
This type of problem occurs frequently in robust control and is known as a generalized 
eigenvalue optimization problem (GEVP) [3]. The GEVP is not an LMI because (5) is 
not jointly linear in p 2 and P. One simple approach to solving the GEVP is to perform 
a bisection search on p 2 , but there are more sophisticated methods available; see for 
example [2]. 


Remark 13. These results may also be carried through in continuous time. In that case, 
an equation analogous to (4) must be satisfied for G(s — A) for all uj £ [0, oo), and can be 
verified by finding P = P T , A > 0 such that 


'A T P + pA- 2A p 

PB 

+ A 

'C T 

b t p 
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P T _ 


M [C 


D] A 0 


5 Examples 

In this section, we show some classes of nonlinearities that can be described by p-IQCs and 
therefore used in Theorem 8 to prove robust exponential stability of an interconnected 
system. In the case where p = 1, these p-IQCs reduce to standard IQCs [13]. This 
class of IQCs will be constructed for SISO systems, but they may be adapted for square 
MIMO systems where the nonlinearity is of the form diag({A}) for a scalar A, with little 
modification. 
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5.1 Pointwise IQCs 


A nonlinearity A satisfies a pointwise IQC with a factorization (*1 f,M) if zjMzk > 0 for 
each k. In other words, the IQC holds pointwise in time. In this case, A also satisfies the 
associated p-IQC for all p < 1. Examples of pointwise IQCs include the 7 norm-bounded 
IQC 


and the [a, /3] sector bounded IQC , given by 


—2a/3 a + f3 
a + /3 -2 


Note that the norm-bounded IQC is a special case of the sector IQC with the sector 
[— 7 , 7 ]. These IQCs hold even if A is time-varying. 


5.2 Zames-Falb IQCs 

A nonlinearity A is slope-restricted on [a,/3] where 0 < a < (3 if the following relation 
holds for all x, y. 

(A(x) - A (y) - a(x - y)) T (A(x) - A(y) - fi(x - y)) < 0 . 

This relation states that the chord joining input-output pairs of A has a slope that is 
bounded between a and /?. This class of functions satisfies the so-called Zames-Falb 
family of IQCs [7, 21]. We give the definition below. 

Proposition 14. A nonlinearity A that is static and slope-restricted on [a,/3] satisfies 
the Zames-Falb IQC 

[-a/3(2 -H-H*) a{l-H)+/3{l-H*)] 

[a(l-H*)+/3(l-H) -(2 -H-H*) \ W 

where H{z) is any proper transfer function with impulse response h := (ho,hi,...) that 

satisfies ||ft||i < 1 and hk > 0 for all k. 

Proof. See for example [7]. ■ 


Remark 15. 


The Zames-Falb IQC ( 6 ) admits 


the factorization 


A! = 


7?(1 - H) 
—a 


— (1 — H) 

1 


and 



1 

0 


In general, for a given fixed p , only a subset of the Zames-Falb IQCs will be p-IQCs. 
We now give a characterization of this subset. 

Theorem 16 (Zames-Falb p-IQC). Suppose A is static and slope-restricted on [a,/?]. 
Then A £ IQC(II( 2 ), p) where II is the Zames-Falb IQC ( 6 ) and H satisfies the additional 
constraint 

OO 

J2p~ 2kh ^ 1 ( ? ) 

k =0 


Proof. The proof involves rewriting the IQC as a discrete-time sum which can be split 
into parts that can separately be shown to be nonnegative. See the Appendix for the full 
proof. ■ 



5.3 Multiple IQCs 


Much like how multiple IQCs can give more precise L 2 gain bounds, multiple p-IQCs can 
give more precise convergence rates. We present numerical examples with both pointwise 
and dynamic p-IQCs. Consider a stable discrete-time LTI system G(z) in feedback with 
the sigmoidal nonlinearity A(x) = 6arctan(x). This interconnection is shown in Fig. 5. 



Figure 5: LTI system G in feedback with the static nonlinearity A(a;) = 6arctan(;r). 


Since this nonlinearity is static, in the [0, 6] sector, and [0, b\ slope-restricted, it satisfies 
the following p-IQCs 


(norm-bounded) (8) 

(sector bounded) (9) 

,2 ki z -k + I (off-by-fc Zames-Falb) (10) 


n n ( 2 ) := 

"6 2 0 

0 -1 

n 0 ( 2 ) := 

'0 6 ' 

b -2 

life ( 2 ) := 

0 

6(1 - p- 


6(1 - p 2k z~ k ) 
k ) -2 + P 

where we can choose any k = 1,2,.... 


A tight bound. For our first example, we analyzed the following LTI system 1 

(* + 1)(10j8 + 9) 


G 1 (z) := - 


(2z - l)(5z- l)(10z- 1) 


(ID 


We solved the feasibility LMI (5) using MATLAB together with the CVX package [5, 6] 
to find the fastest guaranteed rate of convergence. We searched over positive linear 
combinations of subsets of the IQCs (8)-(10). Fig. 6 shows the rate bounds achieved as 
a function of which IQCs were used. For the particular choice 6 = 1, Fig. 7 shows sample 
state trajectories. 

The true exponential rate can be found by linearizing the system about its equilibrium 
point. Namely, A(a;) ss bx. Formally, this is an application of Lyapunov’s indirect 
method [10, Thm. 4.13]. The result is that the decay rate should correspond to the 
maximal pole magnitude of the closed-loop map G(z)/(1 — bG(z)). We display the true 
exponential rate as the dashed black curve in Fig. 6 and Fig. 7. 

For this example, the p-IQC approach yields a tight upper bound to the true exponential 
rate when we use a combination of the sector and off-by-1 IQCs. 


A loose bound. The p-IQC approach does not always achieve tight bounds as in the 
previous example. Consider the same problem as before but this time using 


G 2 (z) 


2z - 1 

10 ( 2 ^ 2-2 + 1 ) 


( 12 ) 


lr This example was inspired by the continuous time example given in [17], which showed that adding 
more IQCs yields better L 2 gain bounds. 
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Nonlinearity parameter b 


Figure 6: Upper bounds on the exponential convergence rate p for the system G\(z) given 
in (11) in feedback as in Fig. 5. A tight bound is achieved using two p-IQCs. 



Figure 7: State decay over time of the system G\{z) in feedback as in Fig. 5 with b = 1 
for various initial conditions x 0 G [—15,15]. The dashed black line is p k , where p = .7058 
is the true rate at b = 1 in Fig. 6. 


The rate bounds for various p-IQCs are shown in Fig. 8. This time, we again observe 
that using more IQCs achieves better rate bounds, but the bound is not tight even after- 
using six IQCs. 


6 Conclusion 

We presented a modification of IQC theory that allows the certification of exponential 
rates. Although we only gave the p-IQC specialization for pointwise and Zanres-Falb 
IQCs, the concept can in principle be extended to other IQCs such as, for example, 
uncertain time delays or slowly varying systems [13]. As the dictionary of p-IQCs is further- 
populated, the applicability of the technique outlined herein would be correspondingly 
expanded. 
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Figure 8 : Upper bounds on the exponential convergence rate p for the system £ 2 ( 2 ) given 
in (12) in feedback as in Fig. 5. As we include more p-IQCs, we can certify tighter bounds. 

A Appendix 

A.l Proof of Proposition 5 

Suppose the interconnection of Fig. 3 is stable. Then there exists some K > 0 such that 
for any choice of the signals e and / and for all T, 

T T 

EOKH 2 + IHI 2 ) ^ ^EOMI 2 + ll/fell 2 ) (is) 

k- 0 k—0 

The proof will follow by carefully choosing e and / to transform Fig. 3 into Fig. 1. To 
this end, note that (A, B ) is controllable by assumption. So there exists a finite sequence 
of inputs uo,..., u n -1 and corresponding outputs ijq, .... y n -i that drives the state of G 
from £ 0 = 0 to £„ = xq. Therefore, if we set 

p~ k u k 0 <k<n (— P~ k 'Vk 0 <k<n 

0 k > n ’ k ~ [O k > n 

then we obtain £„ = xq in the interconnection of Fig. 3. Moreover, fk = e k = 0 for k > n, 
so we may cancel the p blocks using the fact that p_ o p + is the identity operator. It 
follows that for k > n, the two interconnections become identical and therefore = x k - n - 
Substituting into (13), we conclude that 

T n—1 

Edl^ll 2 + INI 2 ) < K E(l|e fe || 2 + ll/fcll 2 ) (14) 

k— 0 k— 0 

The right-hand side of (14) is independent of T, but (14) holds for all T so we must have 
lim ||tUfc|| = 0 and lim ||ufc|| = 0 

k—¥ 00 k—>o o 

For k > n, we have w k = p~ k u k and v k = p~ k y k ■ Therefore there exists some constant 
c > 0 such that 

IKII<cp fc and \\y k \\<cp k 

Now ( A , C ) is observable by assumption, so let L be such that the eigenvalues of A + LC 
are all zero. Rewrite the dynamics of G as 

x k +i = Ax k + Bh k (15) 
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where A := A + LC, B := [LD + B — L\ , and h k := [uj, yJ] T . Iterating (15), we 
obtain 

k -1 

x k = A k x o + A k ~ 1 ~ i Bh i (16) 

»=o 

Since all eigenvalues of A are zero, A is nilpotent and so A n = 0. For k > n, (16) therefore 
becomes 

n —1 

x k = Y An ~ 1 ~ i Bh k - n . 

i=0 

We can now bound the state using the triangle inequality. 


k -1 




... as /r! ^ \\hi 


i—k—r 


<7 c 


p“ n - 1 

1-/3 


and this completes the proof. 


A.2 Proof of Theorem 16 


We will prove this general result by first considering the simpler case where the slope 
restriction is on [a,/3] = [0, oo] and H(z) = p 2j z~ j . Note that this choice trivially 
satisfies (7). In this case, the II from (6) becomes 


n = 


0 1 - p 2j z~ f 

1 - p 2j z~ j 0 


(17) 


where z denotes the complex conjugate of 0 . We call (17) the “off-by-j” Zames-Falb IQC. 
We would like to show that A £ IQC(n(z), p). Appealing to Definition 6 and Remarks 10 
and 15, this amounts to proving that 


Y p 2ku k(vk - P 2 °yk-j) > 0 

k =0 


(18) 


We will prove (18) by borrowing the approach from [12]. If A is multidimensional, we 
require that A be the gradient of a potential function [7]. By the assumption that A is 
slope-restricted on [0, oo], we have 

(A(x) — A(y)) T (x — y) > 0 holds for all x, y 


In other words, A is monotone. Now define the scalar function g such that Vg = A. By 
Kachurovskii’s theorem, g is convex and we have 

g(y) > g(x) + A(x) T (y - x) for all x, y 

Moreover, setting (x,y) >->• (y fc , 0) or (x,y) >->• {y k ,y k -j) leads to the two inequalities: 

ujyk > g{y k ) 

ul(y k - yk-j) > g(yk) - g{yk-j) 
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(19) 

( 20 ) 



We will assume for simplicity that g(x) > 0 for all x. Substituting (19) and (20) into the 
left-hand side of (18), the partial sum from 0 to T is: 

T 

Y P~ 2ku k(Vk ~ P 2j Vk-j) 
k—0 

T 

= Y P~ 2k (( l - P 2j ) u lvk + P 2j ul{y k - Vk-j)) 
k—0 
T 

> Y p ~ 2k (t 1 ~ + p2j (9^ -PiVk-j))) 

k—0 

T 

= Y p ~ 2k ( 9 ^ ~ p 2j 9(yk-j)) 

k—0 

T 

= Y P~ 2k 9(Vk) > 0 

k-T-j+l 

Since each partial sum is nonnegative, the infinite sum (which must converge) is also 
nonnegative, and therefore we have proven (18). For the case of general [a,/3], observe 
that II from (6) may be factored as 


' p -1' 

T 

0 

1 - H(z)*' 

" P -1' 

—a 1 


1 - H(z) 

0 

—a 1 


So the above proof holds verbatim if we make the substitution (y, u) >->• (/3y — u,u — ay). 

Now we consider the case of a more general H(z). Suppose H(z) = hkZ~ k where 

hk satisfies (7). Then, 

OO 

1 - H{z) = 1~Y h * z ~ k 

k =0 

OO OO 

= 1 ~Y P~ 2k bk + Y P~ 2k bk (! - P 2 k z~ k ) 

k =0 k =0 

N_ v ✓ 

C 

OO 

= c(l-H s ) + Yp~ 2k hk(l~H k (z)) , 

fc=o 

where H k (z) = p 2k z~ k and H s = 0. Note that H k (z) corresponds the off-by-fc Zames-Falb 
IQC, which we proved above is a p-IQC. Also, H s corresponds to the sector IQC, which 
is also a p-IQC. Now note that the general Zames-Falb IQC (6) is linear in I — H and 
/ — H*. Therefore, II(z) is a positive linear combination of p-IQCs and must therefore be 
a p-IQC itself. ■ 
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